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I. mTRODUCTION 


The peach crop is a major contributor to the agricultural output of the 
Southeastern United States, particularly the state of Georgia. However, the in- 
tensive culture of peaches over a long period of time leads to the so-called 
"peach decline" . Over the last few decades, the average productive life of peach 
trees in Georgia has declined from approximately twenty years to from five to 

ei^t years. 

Peach decline has been investigated in Georgia since 1929 The premature 
death of peach trees in the Southeast has been attributed to such disorders as 
root rot, peachtree borer, nematodes, bacterial canker, fungi, cold damage, 
and virus diseases. 

In 1972, between 100, 000 and 200, 000 trees died, and hence in that year 
aerial photography of several orchard test sites was obtained. Previously, the 
occurrence of tree decline due to the presence of three diseases, three insects, 
and one mite had been detected by the examination of color infrared aenal photo- 
graphs obtained at elevations of 600-4500 feet. ^ 

Color infrared photographs obtained in 1972 were analyzed using the 

2 

General Electric Image 100 System By displaying on a color television moni- 
tor, categories corresponding to selected density levels on the film are readily 
located. With this method, trees in various stages of decline and also different 
segments of the same tree fall into different categories. 

Consequently, analysis by digital computer of multispectral aerial photo- 
graphy was undertaken. The purposes of this study were to determine the 
feasibilities of orchard inventories (enumeration and sizing of the crowns) 
and of determining the state of decline of trees based on an analysis of their 
multispectral reflectances . 

This report presents the results of the analysis of four band multispectral 
photography by digital computer The major topics presented are (i) digital 



correction of intensity variations, (ii) enumeration of trees by sizes, and (lii) 
classification into healthy and declimng categories. 

The aerial photography was undertaken by NASA/MSFC, and ground'truth 
data was furnished by the USDA Southeastern Fruit and Tree Nut Research 
Station, Byron, Georgia. The aerial imagery was obtained over the USDA test 
sites in Bibb , Crawford and Peach Counties in the state of Georgia The loca- 
tions of these test sites and the USDA Research Station are shown in Figure 1.1. 

O 

The photographs used were made by the MSFC I S four band camera system, 
which exposes four frames of film size 90x90 mm. through filters which trans- 
mit in the blue, green, red and infrared spectral bands. Digitization of the images 
and computer analysis were performed at Marshall Space Flight Center by Com- 
puter Sciences Corporation. 
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Fig 1* i ’Locations of peach orchard test sites and USDA Research Station 
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n. DIGITAL FILTERING TECHNIQUES 


2.1 T.w,o Dimensional Linear Filtering 

The first topic to be discussed is digital correction of intensity variations, 
by means of digital filters. This method is feasible as a procedure for removing 
variations which have spatial rates of change differing from those of the objects 
of interest in the scene. This is particularly true of an aerial photograph of a 
peach orchard, in which the objects of interest are of a consistent size. 

The input/output relation of a linear digital filter is of the form 
L 

Zf — 23 Sk ^ i+k 

k=-L 

where [x] is the input data sequence and [Z] is the output data sequence. It 
may be seen that an output data value is a sum of the neighboring 2L input values, 
plus the input data point, each multiplied by a filter weight gj^ The implementation 
of a filter consists of selecting the filter weights, followed by evaluating the in- 
put/output relation at each data point. Because of the large amount of data in a 
digitized image, it is essential that the implementations be fast in order to make 
digital processing practically feasible 

A very fast implementation is obtained when the filter weights are given 
the value one, and the evaluation is done recursively In this case, the filter 
output is a moving sum over the input data, and in advancing by one data point, 
it is necessary to add only one new sample and subtract one sample which has 
left the filter length of (2L+1) samples The following diagram illustrates* 


i-L-1 

i-L 



l+L-1 



1-L 

i-L+l 


i+L-1 

i+L^ 


= Zj_i + Xi+L - ^i-L-1 
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Also, in the actual filtering of a digitized image, two dimensional filtering is 
performed This involves extending the summation over all neighboring points 
in a region of dimensions (2Lx+l) (2Ly +1) The output data values will occupy 
the same magnitude range as the input if the average of the (2Lx+l) (2Ly+l) terms 
IS taken. 

When the two-dimensional moving averages are applied, data variations 
which occur over a span of less than 2L+1 data points tend to be averaged out, 
while more slowly changing variations remain Hence, the operation has the 
effect of low pass filtering Denoting the moving average low pass filter oper- 
ation by LP we can write 

Z = LP (K) 

Since the inverse of a low pass filter transfer function is that of a high pass 
filter, a high pass operation may be obtained as 

Z = HP (X) = [1 -LP] (X) = pC) - LP (X) . 

This IS equivalent to subtracting the moving averages (the slowly varying inten- 
sity variations) from each data point constant bias may be added to restore 
the data to its original range) . 

Band pass filters may be constructed from combinations of low and high 
pass. A filter which passes in the midrange is given by 

BP = LP (Li) - LP (L2) 

where Lj^ > L 2 . In this case, LP (L^) removes the extreme high frequency 
variations, while subtracting the output of LP (L 2 ) removes the extreme low 
frequency variation. A filter which attenuates the midrange is given by 

BP = LP(Li) HP (L 2 ). 

In this case, the high and low ranges are passed. 



2 . 2 Digital Filtering of a Digitized Peach Orchard Aenal Photograph 


Initial experiments were performed using a green band image obtained at 
6000-feet altitude. A rectangular orchard consisting of approximately 37x84 
rows of peach trees was digitized to 64 levels spanmng a range of 0 to 2D in 
photographic density. The digitizing spot size was 50 microns, yielding for the 
digitized scene an array of 850x1800 pixels 

Visual observation of the image indicated large variations in scene inten- 
sity, although maintaining distinct contrast between the tree crowns and back- 
ground. The histogram of the scene density is given in Figure 2 1. The histo- 
gram verifies the wide range of densities present in the image, and lack of 
uniformly separable intensity ranges corresponding to tree crowns and back- 
ground 
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Figure 2.1. Histogram ol densities present in the peach 
orcliard image. 
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As a first step in the processing of the image, it was decided to attempt 
the removal of the background intensity variations . The possibility of perform- 
ing this operation arises from the slower rates of spatial change of the back- 
ground intensity variations compared to the variations due to the tree crowns 
This fact IS illustrated in Figure 2 2, which is a trace of the density variations 
along a row of trees as the background density varies widely. The ten intensity 
peaks correspond to ten tree crowns which are dark in appearance on the photo- 
graph The overall intensity variation is evident in the plot, with a region of 
low intensity near pixel 100. 

In this case the variations of intensity do not follow a positional dependence 
which would allow fitting to a function of position. Hence the moving averages 
(low pass output) must be computed for each point in the digitized image. If 
the size of the moving average region is chosen correctly, the low pass output 
will be a measure of the low frequency variations in image density. This out- 

t 

put at each point in the image is subtracted from the original data value and a 
bias IS added to restore the overall average density of the output image. As 
indicated previously, the result of subtracting a low pass filter output from the 
original data is characteristic of high pass filtering. Figure 2.3 is a trace 
along the row of pixels shown in Figure 2 2, after the image has been high 
pass filtered Twenty-five neighboring points were taken as the filter length 
in each direction, yielding a total area of 51x51, including the data point itself 
It can be seen that the high pass output retains the detailed structure of the 
original data while significantly reducing the slow drift exhibited in the data. 

A further refinement was attempted using the high pass filtered data 
This data was low pass filtered, using a filter length small compared to the 
dimensions of the tree crowns Thus, in this case, the variations arising from 
the tree crowns are considered as slowly changing and are relatively unaffected, 
while smaller patches of vegetation and nonuniformily within tree crowns are 
reduced. The net result of these operations is a band pass filtered image, in 
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Fig. 2, 2, Density variations along a row of trees in the original data. 
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Fig. 2.3. Density variation along a row of trees after high pass filtering. 
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which only those objects which axe the size of tree crowns retain their original 
intensities. Figure 2.4 shows the original data sample and the low pass and 
band pass outputs . 

Examination of the original data scan shows that a separation of tree 
crowns and background is not possible for any value of threshold. However in 
the filtered output image all of the tree crowns can be separated from their 
surrounding background by a density threshold of approximately 25 (on the 0-63 
scale of scanner output). This is further demonstrated by the histogram of the 
filtered data (Figure 2.5), which has become bi-modal, with a minimum occuring 
at the value 25 . 



Figure 2.5. Histogram of densities in the peach orchard image 
after digital filtering. 
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Original image and result of background removal following 
digital filtering of a peach orchard test site. 
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Using a threshold of 25, a digital image can be generated in which the 
tree crowns remain while the background is set to zero density. However, as 
the densities of the tree crowns have been altered somewhat by the filtering 
operation, a more accurate result is obtained if the filtered image is used as 
a mask to be compared with the original image. In this manner, the background 
in the original image can be set to zero density. A resulting image obtained in 
this manner and the original image are shown in Figure 2.6. The histogram of 
the tree crowns only is shown in Figure 2.7. The symmetry of the distribution 
indicates the exclusion of the bacl^round pixels. 



Fig. 2. 7. Histogram of tree crown densities [ following background removal] 
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ni. DETECTION TECHNIQUES 


3.1 Introduction 

The problem of determining the positions of specified objects in a digitized 
image array may be approached by comparing, in some way, a local region in 
the image with a mathematical template which represents the searched-for objects. 
A match is considered to have occurred when the similarity between the image 
and the template reaches a certain level, in terms of a mathematical criterion . 

The traditional, but time consuming, method is a search for peaks in the cross- 

correlation function at image coordinates f i, j] the sum of the product of each 
template term with its corresponding image term. For a tt:mplate T of size 
mxm, with an image array I, the normalized cross-correlation is given by: 


However, it may bo possible to obtain satisfactory results, with a large decrease 

in computation time, by using simple differences between the template and the 
3 

image. 

3 . 2 Template Matching with Simple Differ en ces 

A test area was chosen from the green band image wliich had been high 
pass filtered to remove background variations. The test area of 100 pixels square 
was relatively small, but was so chosen to allow rapid presentation of results 
on the output line printer. A line printer plot of the area is shown in Figure 3.1. 

The template chosen consisted of a circular region inscribed in a square 
of eleven pixels. The template elements representing the disk and the border 
were assigned values approximating the density levels of the tree crowns and 
background . 
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Figure 3.1. Line printer plot of the digitized peach orchard test site. 
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REPRODUCIBILITY OP THE 


The template was centered over each pixel in the image, and the differences 
between each template element and the corresponding image pixels were summed. 
This sum, obtained at each point in the image, is termed the cumulative error. 

A printer plot of this error for the test region is given in Figure 3.2. As ex- 
pected, the errors decrease steadily as the template becomes superimposed on 
each tree crown. The pixel position for which minimum error occurs will be 
considered to be the match position. 

Plots of the errors obtained along two rows of the image are given in 
Figure 3.3. The upper plot is the error along line twenty-one of the test site, 
in an area between rows of trees. Consequently, the matching errors are large 
and fairly constant. The lower plot is the error along a line of tree crowns, and 
demonstrates the wide variation in error and in particular the rapid change in the 
error near the match positions (the minima in the plotted curve). 

It is apparent from the preceding figures that an error threshold can be 
determined such that image points whose corresponding errors exceed the 
threshold are far from the locations of the tree crowns. This is illustrated by 
Figure 3.4 which is a printer plot of the errors at only those pixel locations 
having errors below the thresliold. Comparison with Figure 3.1 demonstrates 
the correlation of these low error locations with the tree crowns in the image. 

This type of thresholding to indicate the locations of the tree crowns is prefer- 
able to simple density thresholding because it includes a dependence on shape 
and not only density. 

3 . 3 Error Evaluation Sequence s 

The definition of the cujiiulativc error at each pixel location is the sum of 
the absolute values of the difference between corresponding template and image 
pixels. However, in the present case, the computations can be formulated to 
reduce the computation time. This is a result of the circular symmetry of the 
template, the definition of the template in terms of only two values, and the 
proportion of background area to tree crown area in tlie imagery. Several 
methods of computation were examined in a series of experiments. 
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Figure 3. 2 Cumulative errors obtained at each pixel in the test scene. 
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Figure3,3, Cumulative errors along line 21 and line 1 1 of the test scene. 
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Figure 3. 4,Cuniiilative Error after Thresholding 
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(i) Natural Sequence 

This method is a straightforward summation of the absolute values 
of the differences for all the terms of the template taken in a raster scan 
order. The computation time required on the test site was 91 seconds. 

(All times given are computation times only, and do not include time re- 
quired for input/output operations. Machine used was a 7094/1.) 

(ii) Natural Sequence + Total Error Threshold 

Since the only areas of interest are those with errors below the 
threshold, the accumulation of error may be terminated whenever the 
threshold is reached. This procedure reduced the computation time to 78 
seconds. 

(ill) Error Sort Method - A 

The following two methods are based on sorting of the error terms 
such that the template evaluation sequence corresponds to the errors ob- 
tained, from the highest to the lowest In tins case, the template is shifted 
±1 pixel from a match position with itself along rows and columns, the 

errors in the four cases evaluated, and sorted so that the evaluation 

4 

sequence determines errors in order ot their magnitude The resulting 
sequence is given m Figure 3 5, and the computation time required was 
80 seconds. The method is not advantageous for this particular problem, 
probably due to the small percentage of evaluations which are performed near 
match positions 

(iv) Error Sort Method - B 

This IS a more generalized error sort method. The test scene was 
used as a "training” area, to determine the average error for each term 
in the template as the template was placed at each pixel in the scene . The 
average errors at each point in the temjilate are given in Figure 3.6, The 
smaller errors in the corners are at the template positions corresponding 
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Figure 3 5 Evaluation Sequence for Error Sort Method-A 
(Method No iii) 

to the background of the scene, and are smaller since the majority of the 
scene is background. The evaluation sequence is given in Figure 3 7, 
and the computation time required was 71 seconds. This method has the 
disadvantage of being traimng data dependent, and hence similar, but 
circularly symmetric sequences were tested 

(v) Disk - Background Sequence 

As abetter approximation to the empirical sequence of method No. 
(iv), a sequence was evaluated which spirals inward over the region of the 
disk, followed by an outward spiral over the region of the background The 
computation time was 72 seconds. 
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Figure 3. 6» Average template errors throughout the test scene 
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Figure 3. 7. Evaluation sequence for Error Sort Method-B (Method No.iv] 
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(vi) Increasing Radius Sequence 

The sequence in order of increasing radius evaluates the error for 
all of the template terms representing tree crowns before proceeding to 
the terms representing background. Since the majority of the scene is not 
tree crowns, errors on the average accumulate rapidly and this sequence 
yields a rapid computation, requiring 62 seconds . 

(vii) Elimination of Absolute Value Computation 

In general, the simple difference between a template term and an 
image pixel must be the absolute value of the difference to allow for image 
data values above and below the template values However, in the present 
case of a double-valued template, the template values may be chosen to 
insure that all difference terms have the same sign. For example, the 
tree crowns have high density in the transparencies, and the corresponding 
template region is set at the maximum image value. This is in fact 
preferable in the present case, since the optimal match is obtained for the 
densest tree crown and lightest background. Using this procedure, com- 
bined with the increasing radius sequence and error thresholding , the 
computation time drops to 46 seconds. 

(vii) Recursive Evaluation 

The recursive evaluation sequence takes account of the fact that as 
the template is centered on successive pixels of the image, the majority 
of the error terms retain their previous values. Hence it is only necessary 
to consider the new terms appearing at the leading edges of the template 
regions (both the leading terms of the whole template and each disk- 

- t 

background interface), and the terms being dropped at the trailing edges. 

J 

This method is extremely rapid, requiring 18 seconds, and it is a property 
of the method that the complete cumulative error is obtained 
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3,4 Error Threshold Sequences 


The result shown in Figure 3.4, obtained by thresholding the cumulative 
error, indicates the locations and approximate sizes of the tree crowns This 
result IS preferable to simple thresholding, in which all pixels above the threshold 
would be retained, irregardless of the shapes of the regions retained A further 
increase in computation speed is obtained when a constant threshold is replaced 
by a monotonically increasing threshold sequence In this case the threshold may 
be exceeded when a smaller number of error terms has been accumulated. 

The threshold sequence cui*vo should approximate, in its average slope, the 
curve of cumulative error versus number of terms. An example is given in 
Figure 3 8. The error at test points A and B accumulates so rapidly that they 
are discarded as candidates for match after very little computation. Point C, 
which is at or near a match, is not discarded although the error for all terms is 
accumulated. 




I f 

I RfOION WHERI. 
CALCULATIONS 
ARE ADDED 


35 AO ' 50 


Figure 3 8. Redaction of Calculations with use of a Monotonically 
Increasing Threshold Sequence. 
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At the point, of match, the cumulative error is the sum of the noise due to 
the fact that the areas being compared are different. If this image noise has a 
mean \ and zero deviation, the ratio of cumulative error to mean of the noise, 
E/X , is umty at each point in the sum, and the cumulative error is linear with 
the number of terms . However, an increase in the variance of image noise is 
reflected in an increased rate of accumulation of error, and the threshold curve 
should take on higher values . Barnea and SilvermanS have calculated threshold 
sequences for various deviations. A set of such curves is given in Figure 3.9. 
Barnea and Silverman also present calculated curves based on the probability of 
the cumulative error exceeding the threshold being held constant at each term 
in the sum. A set of these curves is given in Figure 3.10. Further analysis of 
this type of computation is given by Ramapriyan.5 

Error sequences were plotted for the peach test site, using an evaluation 
sequence which starts at the center of the circular template and spirals outward 
The cumulative errors at each term in the sum over a template of size 11 x 11 
pixels are plotted as the solid lines in Figure 3.11. The upper plot is for a non- 
match position, while the lower is obtained for the template centered on the 
upper left tree of the test site The area of the template circular, region was 
nominally 100 pixels, but an actual area of 97 pixels was obtained, due to the 
discreteness of pixels. Thus, after 97 terms of the error sum, there is an 
abrupt change in the template value, and this accounts for the slope change in the 
error sum plots near the 100-th term. 

From the cumulative error plots, threshold sequences can be obtained such 
that the error sums at non-match points exceed the threshold terms, while 
remaining below the threshold near match points. For this test, linear threshold 
sequences were chosen, since the error sequences are nearly linear for the 
imtial terms . The threshold sequences shown as dotted lines increase by 8 and 5 
for each term in the sum. It may be seen from the plots that, for non-match 
positions, the initial term of the error sum will exceed the threshold. The linear 
threshold curve passing through the origin does not make allowance for variance 
m the image noise The smoothness of the actual cumulative error plot justifies 
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CRR0R2 ERROR] 




Figure 3.11. Cumulative error at a non-match and a match point. 

The dotted lines are threshold sequences with slopes of 8 and 5. 
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this assumption. Similarly, it may be seen that, near match points, the error 
IS accumulated for all terms of the template, without exceeding the threshold. 
Hence, an appropriate output function is the number of terms in the error sum 
at which the threshold is exceeded, or the total number of terms when the thresh- 
old IS not exceeded. Printer plots of this output are given in Figures 3.12 and 
3.13 for threshold sequence slopes of 8 and 5, respectively. It is apparent that 
the sums are terminated quickly at non-match points, effecting a large reduction 
in computation time. The actual computation times were 11 seconds and 7 7 
seconds, respectively The regions in which the error accumulates for a large 
number of terms approximate the sizes of the trees in Figures 3.12 and 3.13, 
as do the regions where the error is below a constant threshold m Figure 3.4. 
However, the result in Figure 3 . 13 was obtained with a six-fold increase in 
computation speed. 
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Figure 3. 12 -Plot of Error Sequence Terms for Threshold Sequence Slope 
of 8 Computation time 11 seconds. 
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Figure 3.13 Plot ot Error Sequence Terms for Threshold Sequence Slope 
of 5. Computation Time 7.7 Seconds 
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3 . 5 Enumeration, Sizing and Coordinate Determination 


3.5.1 Techmque 

Previous sections have demonstrated that simple differencing of a cir- 
cular template and a digitized image of a peach orchard produces an error 
function with pronounced minima at the locations of tree crowns . The coordinates 
of a mimmum are readily obtained and may be taken as the coordinates of the 
tree center. However, it is necessaiy to repeat the process of finding a new 
minimum value in the vicinity of each tree crown Hence an approximate position 
of match is found initially by considering the gradient between the leading and 
trailing edges of the circular template positions in the data. 



Figure 3.14 Positive and Negative Gradient terms for a Circular 
Template of Nominal Area iOO Pixels, 
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The template gradient is defined as the positive sum of all terms on the 
leading edge of the circular region, plus the negative sum of all terms on the 
trailing edge The positive and negative terms for the horizontal gradient of a 
template with a circular area of 97 pixels are given in Figure 3 . 14 The hori- 
zontal gradients along lines 13 and 23 of the test site are plotted in Figure 3,15. 
The upper plot is along a line centered on a row of trees, while the lower is be- 
tween rows . 

As a consequence of the facts that (i) the leading edge of the template de- 
fines the image points contributing to the positive terms in the gradient, and (li) 
the tree crowns have digitized values of greater magnitude than the background, 
the gradient values are positive as the template approaches the position of a tree 
crown in a raster type scan. The value of the gradient falls rapidly through zero 
as the template passes over the center of a tree crown. By detecting a similar 
passage through zero in the orthogonal direction of scan, it is possible to rapidly 
obtain approximate coordinates of match 

An alternative description of this procedure is obtained by regarding it as 
a recursive calculation of the sum of the image values corresponding to the cir- 
cular area of the template disk. The gradient terms at each pixel location are the 
terms to be added to the preceding sum to obtain the magnitude of the sum at the 
current pixel location. Thus the gradients are the differential of the sums over 
circular regions centered at each pixel location, and as such pass through zero 
at the extremum of the sum. 

This method of locating the approximate positions of match was chosen 
over others (such as total error threshold or moaoh>nic threshold sequence] for 
the following reasons: 

(i) A recursive-type computation is approximately as rapid as 
a threshold typo, 

(li) passage through zero is easily detected. 
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(iii) the effect of bias in the image density levels on the choice of 
threshold values is eliminated, 

(iv) differentials involving a small number of terms at the edges of 
the template disk yielded inaccurate match positions 

After the approximate match position is found by means of gradient 
sequences, a block of area five pixels by five pixels centered on the approximate 
position is analyzed in detail. In this region, the total cumulative error between 
the template and the image is calculated using a recursive method, and the posi- 
tion of minimum error is determined If this error does not exceed a predeter- 
mined threshold, this position is accepted as the coordinates of the center of a 
tree crown In addition, a tree count is readily obtained. 

3.5.2 Test Site Tree Detection 

The peach orcliaid test site was analyzed in this manner. Figure 3.16 
is an intensity plot of the cumulative errors which were computed over the five 
pixel by five pixel areas centered on the approximate match locations. The in- 
tense positions withm regions of small error are the fmal tree coordinate 
positions, as determined by the minimum total cummulative error. The total 
error is above the threshold for those five by five squares which are intense at 
all positions, as in the leftmost three in the second row At these positions, the 
gradient procedure was sufficiently sensitive to detect the presence of an object, 
but the size or contrast did not match the template very closely, and hence the 
cumulative error was large. 

The test site was anah'zod using total error thresholds of 1000 and 1500. 
The computer printouts of tree coordinates and enumerations are given in Figures 
3.17 and 3 18 

In the first case, the average total error for the 21 matches was 634, or 
5 2 per difference term between the template and the image In the second case 
higher error positions wore accepted as matches, and the corresponding figures 
are 768 or 6.3. Line prmter plots of the test sites with the outlines of the 
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Figure 3. 16. Cumulative errors for regions surrounding the approximate 
match positions. 
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TREE COORDINATES 
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Figure-3r4^7 — T*"*® Coordinates and Enumeration using a Threshold 
of 1000. 
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TREE COOkUINATES 

LINE NUMBER SAMr>LE NUMBER AND TOTAL ERROR 


7 

1C4 

- 618 

126 

- 014 

10 

73 

- 639 



1 1 

31 

- 515 

54 

- 673 

12 

12 

- 547 



2R 

123 

- 33 7 



.30 

73 

- 758 

103 

- 661 

3^ 

10 

-1318 



46 

125 

- 7B6 



49 

104 

- 575 



50 

72 

- 557 



51 

31 

- 571 



5? 

52 

- 636 



6 9 

124 

- 664 



70 

102 

- 602 



71 

72 

- 601 ■ 



73 

53 

- 737 



74 

30 

-1^70 



73 

11 

-12H4 



91 

125 

-1462 



92 

71 

- 4.3 

103 

- 462 

9 A 

5? 

-11A5 



95 

30 

- 603 




IREE COUNT 
NUMBER AREA 

4; 4c 4r 4c 

26 97 


F.LAPSEi' MMH = 12.30 SECONDS 


Figure 3. 18 


Tree Coordinates and Enumeration using a Threshold 
of 1500. 
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detected trees superimposed are given in Figures 3.19 and 3.20. It may be 
noticed that the higher error threshold results in the inclusion of additional trees 
which have weak foliage and exhibit less contrast with the background. The 
apparently paradoxical reduction in computation time as additional trees are de- 
tected is due to the fact that an area occupied by a previously detected tree is 
not searched in succeeding scans, but rather bypassed. 

3.5.3 Programming Considerations 

This section presents a description of the algorithm in greater detail, 
with particular emphasis on methods of computer programming the algorithm. 

In order to rapidly evaluate expressions involving the image pixels cor- 
responding to a template of size nxn pixels, it is necessary to hold n records of 
the image data in core. If the template is to move to adjacent lines of the digi- 
tized image, a corresponding number of additional records should be held in 
core. In the present case the template is allowed to move over five scan lines of 
the image, and hence two additional records of data before and after the n records 
corresponding to the template size are stored. 

Now, the simplest method of referring to the various template terms is 
by numbering them. If the elements of the template are’ numbered along its 
columns (consistent with the storage in memory of a two-dimensional array de- 
fined in the FORTRAN language), the clement (i, j) is located-at location n(j-l)+i. 
The corresponding element of the image data is obtained simply by adding a bias 
term due to the position of the template along the scan line. Thus, one sequence 

of location numbers is determined for each evaluation required, such as the 
gradient terms, cumulative error scans, recursive evaluation sequences. 

Three types of evaluation sequences are used. These are: 

(i) differential sequences used for recursive evaluation of five 
sizes of disk and boixler along rows and columns, 
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Figure 3,20 Detected Trees using a Threshold of 1500, 
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(ii) a sequence over the template in order of increasing radii, used 
for the initial evaluation of the five disk size cumulative sums, 

(iii) a sequence used to flag the area surrounding a match position 
to suppress re-searching that area. 

The initial search for approximate match positions, using the differential 
of the disk sum, is done using the smallest and the largest disk areas, to ensure 
detection of all objects of the appropriate size. (If a large mismatch of sizes 
occurs, the differential sequence does not pass rapidly through zero). A detailed 
analysis is made if the gradient sequences along both rows and columns passing 
through a given pixel consist of two positive terms followed by two negative terms, 
indicating passage through zero. When this occurs, the match coordinates are 
determined as the best position found in a five pixel by five pixel area surrounding 
the approximate match position. 

The cumulative error which is computed is effectively the difference be- 
tween the sum of the image values corresponding to the disk area and those of a 
border area. This is a measurement of the contrast between the circular region 
and its background. In order to compute this measure over five disk areas with- 
out summing over the same area more than once, partial sums are computed over 
annular rings, whose area is the increment in the disk areas. Then sums of 
terms over disk regions and border regions of the five specified areas are obtained 
by summing a small number of the partial sums. 

The additional difference sums for the remainder of the five pixel square 
area are determined recursively, using the sequences of array elements for the 
new and deleted terms as the templates are moved along rows and columns of 
the digitized image. If at the position of maximum disk-border contrast a specified 
threshold is exceeded, the coordinates of the match position and the template 
number (1-5) are output on the printer and sequential mass storage. 

When advancing to the next record of data, it is necessary to rotate the 
records of data held in a two-dimensional array in core by one row, and then 
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read the next record into the last row of the core array. (An alternate method, 
in which the evaluation sequences rather than the imagery data records are 
rotated, was tested. The execution time of this program was greater by one-third). 
Similarly, since the match can occur in any of five lines of data, the match para- 
meters must be placed in arrays of five rows, which must be rotated for each 
record read. 

In order to prevent multiple detection of the same match position, and to 
accelerate the detection process, a specified area surrounding a match position 
is not searched again. This is accomplished by marking the data around the 
match in some manner. This region may be rapidly distinguished from the data 
by the use of negative numbers . If the area to be skipped is large enough that no 
part of the template will fall on this area, the magnitudes of the elements may be 
changed from the data values in order to convey useful information, such as the 
number of pixels to advance along the scan line. This method was used in the 
100x100 pixel test area, where the overall computation time was seen to decrease 
as a greater number of trees were detected (and hence more area was bypassed 
in the search). In general, positions of the template near its edges may fall on 
the marked area, and so the data values are simply negated, so that their mag- 
nitudes remain available. With this method, only the midpoint of the template is 
tested for falling on the marked area, and hence this test is very rapid. However, 
this method limits the dimensions of the bypassed area to the size of the template. 

If a larger bypassed area is desired, the comers of the template must be tested 
for falling on the marked area. 

An additional programming consideration which accelerates the execution 
is the use of a look-ahead test, which cancels the search for passage of the 
differential through zero if there are not sufficient pixels remaining before collision 
with a previously detected match. 
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3.5.4 Large Area Test Site 

The test site known as Bateman's orchard in Bibb County, Georgia was 
photographed by the MSFC l2s camera at an altitude of 3000 feet. The blue band 
image was scanned and digitized to 64 levels in the density range 0 to 2, yielding 
a digitized image size of 600 lines X 850 samples (510, 000 pixels). The orchard 
contained 28X 38 rows of trees, and hence 1064 tree positions. The majority 
of the trees in the orchard have been examined from the ground, and it was 
found that 50.5 percent of the original trees remained alive. Extrapolation to 
the total area of 1064 grid positions indicates 537 live trees present. 

The tree detection algorithm was applied to this digitized image, searching 
for five sizes of trees whose areas ranged from 100 pixels to 300 pixels in incre- 
ments of 50. The corresponding diameters ranged from 11 pixels to 19 pixels 
in increments of 2. The threshold used was 7.0. (A threshold value of 7.0 means 
that, on the 64 level scale of densities, the average density over the disk region 
must exceed that over a ring-shaped border region of comparable area by 7.0), 

In this test site, the tree crowns overlap in many cases, and an important 
parameter is the allowable separation between match positions, which strongly 
affects the execution time due to the bypassing of a specified area surrounding 
previously detected matches . A set of experiments was performed in which the 
bypassed area varied from 9 pixels on a side to 25 pixels. The results are 
presented graphically in Figure 3.21. The computation times do not include I/O 
operations. The machine used was an IBM 360/65. It may be noticed that by- 
passing larger areas decreases the execution time more rapidly than the number 
of detected trees decreases . 

Sample printouts of the coordinates and parameters of each match 
position, and of the summary of match sizes and enumeration, are given in 
Figure 3.22. The five sizes of detected trees, for two values of the area by- 
passed around previous matches, are indicated on the blue band photograph in 
Figure 3.23. 
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Figure 3.21 Computation time and tree count as a function of bypassed -area. 
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Figure 3.22 Sample printout of detection algorithm results. 
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BYPASS SIZE = 2S 


Figure 3.23. Detected trees for two values of the bypassed area 
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3.5.5 Performance Summary 

The performance of the algorithm will be discussed with regard to 
computation speed and detection accuracy. 

As stated previously, the computation speed increases with the number 
of detected trees, due to the bypassing of areas surrounding matches. The present 
scene is not advantageous in this regard, since one-half of the trees are missing. 
Also, the detection rate will depend on the resolution at which the imageiy is 
digitized. It is apparent that, for good detection and sizing, the tree diameters 
should be on the order of several pixels, ranging from eleven to nineteen in the 
present case. With an IBM 360/65, the digitized image was analyzed at an 
average rate of approximately 1350 pixels per second. (377 seconds were re- 
quired with bypass size 25.) Expressed in a different fashion, the rate of analy- 
sis is 3.8 tree cells per second, where a tree cell is a square region of side 
the average distance between rows of trees, which is 18.8 pixels. 

Using the lower image of Figure 3.23, it was determined manually that, 
within the 28 by 38 row area of the orchard, fifteen trees were not detected. 
Subtracting the nineteen trees detected outside of the orchard from the total of 
534 detected trees, we obtain a detection accuracy of 515/530 or 97 percent. 


47 


rV. ANALYSIS OF THE INFRARED BAND PHOTOGRAPH 


It has been known for some time that healthy vegetation is strongly reflective 
in the near infrared spectral band, while vegetation vmder stress from any source 
is markedly less reflective. This trait has been widely exploited in the determi- 
nation of disease and insect infected areas of deciduous forests and citrus or- 
chards. Encouraging results in peach orchard photography have been obtained 

by Georgia fiistitute of Technology in cooperation with the USDA Southeastern 

2 

Fruit and Tree Nut Research Station. In this work, color infrared transparencies 
were viewed throu^ color separation filters by a television camera. The image 
was then analyzed by computer, and it was found that orchard sections containing 
imhealthy trees were discriminated by their infrared reflectances. 

Bateman's orchard in Bibb County, Georgia, was the site of a 9-year experiment 
in various treatment methods. An aerial photograph, with the treatment areas 
marked, is given in Figure 4-1. Herbicides were used to control weeds in strips 
lettered "H". Weeds were controlled by disking in strips lettered "D". Blocks 
of live trees in the disked strips indicate where Fumazone was injected for 
nematode control. The almost perfect stands within the "H" strips indicate 
where both herbicides and Fumazone were applied. 

The infrared band image taken by the MSFC I^S camera was scanned and digitized 
to 64 levels in the density range 0 to 2. In order to establish the existence of 
variations in infrared reflectance as a function of tree vigor, the digitized image 
was thresholded at several closely spaced density levels. The digitized image 
was then rewritten on a film writer with the regions having densities less than 
the thresholds being blacked out. The results for four thresholds are given in 
Figure 4-2. The healthiest trees, having the hipest reflectance, will be blacked 
out at the lowest thresholds. It can be seen that only the healthiest trees have 
densities of 25 or less. Setting the threshold at 28 includes nearly every tree 
within the perfect stands treated by herbicides and Fumazone, while excluding 
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Figure 4-1. Bateman's Orchard Showing Treated Areas 
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Figure 4-2. Infrared Photograph of Bateman's Orchard Showing 
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Figure 4-3. Histograms of Declining and Healthy Trees 
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V. ANALYSIS OF MULTIBAND AERIAL PHOTOGRAPHY 


5.1 INTRODUCTION 

The photography used in this stucfy was obtained with an International Imaging 
2 

Systems (I S) Mark I four-band multispectral camera, using Kodak type 2424 
black and white infrared fQm and type 2420 duplicatmg film. The four spectral 
bands were defined by filters, whose peak transmissions occur at wavelengths 
0.430, 0.530, 0.635, and 0. 800 microns. The transmittance curves are shown 
in Figure 5-1, 

5.2 IMAGE REGISTRATION 

The four photographic transparencies were digitized separately and, hence, 
require digital registration of the Imagery. This was accomplished by choosmg 
five template areas m one photograph and searchmg for match positions in the 

3 

other three images by means of a sequential similarity detection algorithm. 

An error threshold sequence was calculated based on a 10 percent probability of 
exceedmg the threshold at the n-th term of the error summation (see Figure 3-10). 
The templates from the first image were then used to search for matches m each 
of the three remaining images and the match positions were taken as those with 

4 

minimum cumulative error. The four digitized files of data were then trun- 
cated by the appropriate numbers of pixels to produce registered files. These 
are then merged so that the four density levels from correspondmg locations of 
the imagery form a four-dimensional feature vector. These steps are illustrated 
in Figure 5-2. 

5. 3 TRAINING SAMPLE SELECTION 

The data values for all four spectral bands were then selected from the template 
regions correspondmg to 50 declinmg and 30 healthy trees. The total numbers 
of samples obtained were 9906 for the declining and 7718 for the healthy. This 
is a prohibitive number of samples when usmg a nonparametric method of 
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MULTI BAND IMAGE SET TO BE SEPARATED SEPARATED IMAGES TO BE VISUALLY 


CHECKED FOR REGISTRATION AND PREPARED FOR MOUNTING ON THE MICRODENSITOMETER 
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Figure 5-2, Summary of Preprocessmg Requirements 
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determining a set of classification functions. The number of training samples 
was reduced by two methods. 

The simplest method consists of sampling the data values at equal mtervals 
such that a tractable number (e. g. , 500) of samples is obtained. A second 
method consists of obtainmg the four-dimensional histogram of the data in which 
each cell of the histogram represents a unique set of the four data values forming 
a feature vector. The result is 8655 cells for the declining tree data and 6856 
cells for the healihy. A reduction in the number of cells to 993 (healthy) and 
737 (declining) was then accomplished by eliminating all cells having a frequency 
of occurrence of one. This procedure is physically acceptable as these cells are 
at the extreme edges of the distributions and probably represent mixtures with 
soil data values, since the irregular shapes of the tree crowns permit such a 
mixture within the circular template regions. 

Succeeding calculations of the classification functions were coded to use the 
four data components and the frequencies of each cell and thus avoid the repetilaon 
of calculations for feature vectors which are identical. These two methods will 
be referred to as the samplmg method and the cell method m the following 
sections. 

5.4 SUPERVISED CLASSIFICATIONS 

The two supervised classification mefliods which were applied are maximum 
likelihood and sequential linear. 

The maximum likelihood algorithm is based on the assumption that the samples 
within a given class are distributed accordmg to a normal (Gaussian) multi- 
variate density function. 

-(X-m)^ K~^(X-M) 
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where M is the mean vector, K is the covariance matrix, and D is the 
determinant of K. 

Such a distribution is completely specified in terms of a mean vector and a 
covariance matrix for each class, which are determined from the trammg 
samples, hi the classification phase, the probability of an unknown feature 
vector belonging to each class is calculated, and assignment is made to the 
class for which the probability is the greatest. 

Linear diserimmant functions are also determined for each class, using the 
trainmg samples. If the linear discriminant function correspondmg to a certain 
class is evaluated by substituting an unknown feature vector, assignment is 
made to that class when the result is positive. If the assignment is not made, 
the data sample belongs to one of the remaining classes, and successive dis- 
criminant functions are evaluated. This effectively means that the multiclass 
problem is treated as a series of two-class problems, in both the training phase 
(determination of linear discriminant ftmctions) and the classitication phase 
(assignment of data samples to classes). Visualizing the data from different 
classes occupying different regions in the spectral measurement space, one 
can recognize that the outer clusters can be separated from the mner clusters 
and that this process can be executed sequentially. Thus it is necessary to 
order the classes of data such that each can be separated from the remaining 
classes by a linear function (a hyperplane in the spectral measurement space). 
This IS accomplished by computing for each spectral band the totals of the 
separations between all pairs of points in different classes (interclass) and 
within each class (intraclass). The optimum separation is obtained when the 
interclass distance is maximized and the intraclass distance is minimized. The 
coefficients of the discriminant functions may then be determined by maximizing 
a criterion such as the distance of the training samples from the discriminant 

7 

hyperplane. The processing flow is shown m Figure 5-3. Applications of 

8,9 

these algorithms have been discussed in greater detail m other publications. 
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Figure 5-3 • Sequential Linear Classifier System 












Training samples derived by the sampling and the cells methods were used in 
maximimi likelihood and sequential linear classifiers. The digitized scene, 
consisting of 510, 000 pixels, was classified on a pixel -by -pixel basis into three 
classes — healthy trees, declining trees, and soil. These results were then 
further analyzed to obtain the classification of each tree. The tree classes 
were taken as the majority class within the area of each tree as determined by 
the match coordinates and the five sizes of templates. 

The two classification maps are shown in Figure 5-4, obtained using the cells 
method of training samples selection and data from the four spectral bands. The 
classifications of the previously detected trees are shown in Figure 5-5. An 
example of the ou%»ut of the tree classification algorithm is given in Figure 5-6. 
The classification results are given in Table 5-1. When using the infrared band 
photographs only, the two classification algorithms give similar results. The 
density levels assigned to the three classes are as follows: 

LINEAR: 0-30 healthy 31-38 declining 39-63 soil 

GAUSSIAN: 0-30 healthy 31-37 declining 38-63 soil 

5.5 PERFORMANCE SUMMARY 

The goal of the multiband classification procedure is to determine those trees 
which are remaining healthy and those which are in decline. One problem in 
this procedure is the variation within a single tree crown. However, if the 
classification is sufficiently accurate, this may indicate damage in some limbs of 
the tree. An attempt was made to overcome this problem by classifying trees 
as healthy or declining according to the majority of the pixels. It may be seen 
from Table 5-1 that this procedure resvilts in a classification of the 80 ground 
survey trees with an accuracy of approximately 70 percent. Additional ground 
survey information included the statistics of 989 trees, but not the status of 
individual trees. Of 499 live trees, the ground survey indicated 46. 69 percent 
as healthy and 53.31 percent as declining. The computer classifications also 
indicate a majority of trees in the declining category. 
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SEQUENTIAL LINEAR CLASSIFICATION 



MAXIMUM LIKELIHOOD CLASSIFICATION 


Figure 5-4. Classification maps obtained by two methods. 
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Figure 5-5. Classification of the detected trees. 
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Figure 5-6. Output of the tree classification algorithm. 
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Table 5-1. Classification Results 


Classification 

Method 

Training Data 

Training Data 

Orchard Classifications (Percentages) 

Classification 

Accuracy 

Accuracy 

By Pixels 

By Trees 

Time 

By Pixels 

By Trees 

Healthy 

Declining 

Soil 

Healthy 

Declining 

Soil 

(Seconds) 

Sampling Method 










Maximum 

Lilcelihood 

73.7 

76.25 

18.62 

31.14 

50.23 

37.3 9 

62.43 

0. 18 

564 

Linear 

70.9 

67.5 

14.88 

33.29 

51.83 

47.01 

52.27 

0. 73 

126 

Cells Method 










r^laximum 

Likelihood 

79.1 

77.5 

16.48 

43.04 

40.48 

47.00 

53.00 

0.0 

524 

Linear 

76.8 

67.5 

11.24 

26.73 

62.02 

27.15 

71.35 

1.50 

71 

Difrared 

Maximum 

Likelihood 

62.3 

68.75 

17.09 

32.40 

50.51 

41.01 

42.13 

16.85 

183 

Infrared Linear 

63.2 

80.0 

17.09 

39.17 

43.73 

39.70 

51.12 

9.18 

47 

Average 

71.0 

72.9 

15.90 

34.30 

49.80 

39.88 

55.38 

4.74 


Ground Survey 






46.69 

53.31 




The infrared band classifications were performed since it is known that decline 
would be indicated primarily in the infrared spectral range. However, the 
accuracy is diminished in this case. This is probably an indication of the fact 
that the spectral band being used is in the very near infrared. It is anticipated 
that better results could be obtained with better infrared response in the detector 
used. 
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APPENDIX A 


A FORTRAN listing of the complete detection package follows. The core 
required on the IBM 360/65, when using standard FORTRAN input and output 
tapes, is 148 K bytes [ K= 1024] . 
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c. 

c 

X. 

c 

X. 

c 

X 

c 


c 


01HEK5ICN S(2V,350». THPLT(2S,25» 

inf.trAi TMPiT 

DATA hR.IAi^EAl ,MbYP, THRESH, N i^cCSfNSAMP / 2 5 , 53 , 7 , 7. 0 ,60C, 850/ 


NR - LENGTH OF A TEMPLATE SIDE 

lAREtl - TEMPLATE AREA INCRFMEA^J 

MBYP - LENGTH OF A bYPASSEG AREA SIDE 

ThftF<;H. - THQP«:H ri n m AvFati.f r>T<:K and BHRnFR DENSITY DlFFFRFNfF 

NPECS - NUMBER OF RECORDS TO BE READ 

NSAm P - NUM.atJL.XF SAMPLES PER RFfORO 

rail ncF cMr is, tmpil,_nr, t abf ai M aY..E.l 

RFl.INf; If) 

REWIND 11 

rail nETPCT Ni^, lABFai. thrfsh. NRFfS . NlSAMPl 

END FILE 11 

iIH£ 

END 
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c 

5L6ROUTINE OSEQNC (ST, TMPLT. NR. lAREAl. M0YP» 


c 

c 

DEFINE EVALUATICN SEQUENCES 


c 

INTEGER ST(1 ) . HSN( 25.10) . VSN(25.10). HSD(25.1C). VSD(25,10> 


OIMENSIC.VJ NSClOIr ISEQ(625 ). KAREAdll, AS1ZE15.2). HSEO(SOO) 
LPCICAL T-^PLKNR.MRI 

r 

COHHON /SEwNCE/ NS ,HSN , VSN ,HSO , VSD , I SE U, K Ak E A, AS 1 ZE , MS , MS EU 


c 

c 

FIND ELEMENTS PRESENTLY DISK, PREVIOUSLY bACKGRGUNO, AND VICE VERSA 
ASSUME A SEARCH AREA OF SIZE (NR«-4) LINES X NR SAMPLES 


c 

COMPUTE VERTICAL MOTION ELEMENTS BY TRANSFORMATION I » J, J = NR-dl 


N « C 

nn 4-? M = i .? 

DO 43 N2>1 ,5 
N » N ♦ 1 

N3 » G 

M4 « 0 

c 

c 

OBTAIN TEMPLATE M 1 TH DISK AREA MAREA* 


c 

lAREA « NI * (N2^l) • lAREAl 


c_ 

CALL TEMPLT (TMPLT, lAREA, NR» 


DC 32 ! 

DC 32 J»3.SK ! 

IF (J.EC.OI GO TO 44 j 

IF (JaEC.NRi GQ TO 33 


IF ( .NOT. TMPLTU ,J). AND. TMPLT( I,J^1) » GU TO 35 
IF 1 .MCT. TMPLT 11 .J»H. AND. TMPLTI I . JH GQ TO 3fe 

i 

1 

GO TC 32 * 

44 IP (Tf*PLT( I ,J*1) ) GC TO 35 

GO TO 22 ; 

33 IF (TMPLT( I.JI 1 GO TO 36 ! 

C 

GO TO 32 

• 

c 

c 

COMPUTE NEW DISK ELEMENTS 



36 N3 « N3 ♦ 1 

HSMN3.NI « INR*4)*|J-1) ♦ I ♦ 2 

• 

VSMN3iNI » 4 Nk«-4) *(NR-I 1 ♦ J ♦ 2 
GC TD 32 

c 

c 

COMPUTE DELtTED DISK ELEMENTS 


c 

35 N4 * N4 ♦ 1 


HSC(K4,tgi » (NR*4I*(J-1I ♦ I ♦ 2 
VSr(K4.N» . (NR*4)*4NR-I) ♦ J ♦ 2 

32 CONTINLE 
43 NS(M « N3 

c 

c 

DEFINE ISEQ IN ORDER CF RADII BEGINNING AT THE DISK CENTER 


c 

J1 . 0 

. 


no ?•> .i,»i ,f'H 

. . . . * — 
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Jl « J 1 ♦ 1 

ST( Jl) » II **2 ♦ JJ**2 

?? TSFQl.ni - «l .1^-1 1 t T« » ? 

CALL SCRTSL (ST, ISEw, I, NR*NR» 

C 

C COMPUTE SEQUENCE USED FOR DEFINING REGION SURROUNDING A MATCH 

C 

MID * (NR+n/2 

MPAM . HrtYP / 2 

IRECl * “ID - 2 

1R FC2 « M3 ♦ MKAD ♦ 2 

ISAKPl * HID - MR40 

1 SA MP2 « MID ♦ MRAD » 1 

MS » P 

DO IP IW»1REC1 .1REC2 

DC 10 J»»*l SAKPl , I SAPP2 

MS » MS ♦ 1 

10 MSEQ(MS) * (NR*A)*(JW-U 1 W ♦ 2 

C 

C COMPUTE DISK AND BORDER SEQUENCES AND AREAS 

I! C 

DC to 1=1,11 

K/ikFfll M » MtlMlAREAl 

60 CALL TEMPLT (TMPLT, ^AREAm, NR) 

DO 5 NS12E»1.S 

ASUEINSIZE.lt = KAREA(NSIZe) 

S AS1ZF(NS17P.2) « K AR F A ( 2 »N S I ZF ♦ 1 ) - KAREAiNSIIEJ 

RETURN 

EM) 
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c 

SUBRCUTINE TEMPLT (TMPLT. lAREA. NHt 

LOGICAL THPLT(NR.NR) 

£ 

K2 « lAREA / 3.1A15S265 

lAREA » n 

DO IICC 

At « 1 - <NR»1 I/? 

DO IICO J»l,NR 

AJ » J - (NR»1>/? 

TMPLH I ,Jt « .FALSE. 

IF (A1»»2»AJ»*2.GT.R?) GO TO 1100 

lAREA > lAREA ♦ 1 

TmplTU.JI « .TRUE. 

non CONTINUE 

RETURN 

END 



c 

<;nflaniiTiNF nPTFfT is. st, mr, mri , iarfai^ thrfsh, waFfs^ wsampi _ _ 

C 

r nPTFrT ripruiaft KFr-irKiS r.iMTKAST rPMPahisnM 

C 

IWTFf.FR SI ,MR1 ,NS awpi . STIll, SAMPl5,2ni. Sl/£l5.2aj_. DELHt DELVt.. 

.HSN(25,n», bSD(25,10), V5NI25,10», VS0(25,10» 
fMMFIJSirrj MSHIAl. NTESTIsl. NSFuISl. NMATI51. n I F F ( 5 . 20 1 . NSI2.E 

.151, NSm», ISEU 625), KAREA(ll), ASIZE(5,2», MSEU(500» 
n«TA MFST y?*1 , fJSFT /Of 0, 2. 1/ 

CCKKON /SEUNCE/ US ,HSN , V SU f H SD . VS 0 , I SE U t KAK E A t AS UL irtS • HS Ed 

c 

IPO FORMAT ( !H1 ,20x, 'COORDINATES OF MATCH POS 111 ONS * /2 1 X , 30 ( ' • • » //5X , 
.'LIUF SA*'.PLE NUMBER. KFL4TIVE AREA. CONTRAST'//! 

101 FORMAT ( I 1 D,3Xf 8 (I 6 ,13 ,F6. 2 M 

102 FORMAT I I i X , E I t U , t 1 , FA. 2 1 t 

103 FORMAT 1 IHl ,30X , 'COUNT OF MATCH POS ITI ONS ' / 31 X , 241 ">• ) // 2 1 X, 'N UMB E 
.R NTMIMAl ARFA ♦ DIAMFTFR ACTUAL P1X=L AREA ♦ D 1 AM ET Ek* / 21 X. 

.6('*'l,3x,12('*'!,3X,8|'*'),3X,17('*'),3X,d('*'ll 

104 FORMAT 1 / 1 2A . I 12 .F 1 3 . 1 . 2 I 1 5 ! 

c 

PR iNT 1 rn 

MIDRCW * NR/2 ♦ 1 

MID * ( NR»4» *( Ml CROR-1 1 ♦ MIDKOW ♦ 2 

NREC2 * NRECS - NRl ♦ MIDROW - 1 
NSAMP2 » NSAM? - MIDROW 

C 

r INITIALIZE SIZF COUNTER. MATCHES PER RECORD COUNTER. AND ARRAY S 

c 

f)1 2 T *l.«5 

NSIZFC I ) = 0 
2 N M A T I I » * 0 

DO E JW=1,NK1 

S RFAlJ (10! ISIJW.IFLl. IEL*1.NSAMPI 

c 

nn 1 

C 

IFL = MrjR'JH 

JS * 0 
f]5FQ( n » ^ 

MSEC(5) * 0 

C 

5330 CCUTINLE 

lEL = lEL ♦ 1 

IF ( IEL.0T.NSAMP2! (.0 TO 1111 

IS r IS ♦ NRl - 

C 

IF I ST1JS»MID!.GF.0I GO TC 80 

N5EU( 1 ) = 1 
M <; F c 1 s 1 = n 

GG TO 5330 
C 

C CANCEL SEARCH IF INSUFFICIENT PIXELS AHEAD TO SATISFY DIFFERENCE SEO. 

f . _ 

80 NPIX » 3 - MAXO (NSEOm, NSEQ(5II 

JSl « JS > NP1X»KR1 


IF MTI r.n to 
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c 

c 

c 

SEARCH FOR »*A<IHUM IN SOM OVER CIRCULAR REGION 3Y SEARCHING FOR 
n IFFFRFNCE SEJIJFNfF WITH SIGNS » ♦ - - 


1 

c 

ft A 

on AO Ks1 fA 


1 



M a f4S(NI 
nFLH X n 




A5 

00 45 JSEQ*1,N3 

IlFIH r CFIH ■» T t B < t S T t .1 <;♦ H <1N ( J S F Q . N 1 H - 1 AR S 1 ST IJ S HI J S FQ .N >11 



NSEiJ(N) s WSEJ(N> ♦ 1 

IF fDELH»NTEST(NSEaiM 1. LT.Q» GOTO 41 

IF (NSEC(M.EQ.4> GO TO 1250 

GO 


A1 

AO 

NSECIN ) * NSET(NSEQIN) 1 
CONTINUE 



c 


GO TO 5330 



c 

c 

COMPUTE DIFFERENCE SEQUENCE ALONG COLUMNS 




1250 

NSEQ( 1 ) a 0 

n<;fo(si : n . 



00 2100 N*lf5t4 

JSl « JS - 2*MR1 - 1 

N3 * NS(^» 

CO 2C10 K»l.i 

DELV « 0 

DO 3S JSgQ^l«N3 


35 

OELW - DELV ♦ IA6S(ST( JSl^’VSM JSEUtN))t-IABS(ST(JSl^VSC(JSEQ.Nm 
IF IDELV'-NTESTIM. LT.O) CO TO 2100 




2C10 

JSl s JSl ♦ 1 

CALL MATCH 1ST. NRl . N MA T . SAMP. SI2E. DIFF. THRESH. lEL. JSl 




21CQ 

GO TO 5330 
CONTINUE 




iin 

GO TO 5330 

CCAJINLE 



C 

c 

RnTATP ARRA/ £ BY 0 ^P ROW 


1 

c 


nn #1 Jux?.M;ii 


* 


L. 

DO e IEL*1.NSAMP 

S( Jh-1 .!EL» « SI JW.IELI 



c 

c 

OUTPUT L 1 5T 5 QF MATCH COD RDl N AT £S > READ MEXT DATA RECORD. 



c 


_Jf IhMATt li.FO.QI C.n TO PCCQ 



NH * NmaT( 1 ) 

sN_^AICH_4 AH 


2121 

on 2121 1 = 1 , NM 

NS I7F ( <; T ZE n « 1 ) ) a N5I ZE ( SIZE ( 1 . 1 ) > » 1 



IP (MW.GT.3) GG TO 800 

PR 1 KT 101 ^ NR£C « ( SAMP ( 1 « 1 1 • S IZ£ (A ■ 1) . l» 1»NH] 


fl-CQ_ 

GO TO 85 

PRTKT IPl. NPFC. ISAMPfl.n. SlZEIl.il. DIFFIl.ll. I«1.8I 




55. 

PRINT 102. (SAMPIl.I), SUEd.I), OIFFd.II. I«9.NM) 
_hFITE dll NFEC- NM. (SAMPll.d, SI2Ed.II. I»1.NMI 




20C0 

CCNTINLE 

-RFAn I ICl ISINRl.lEU. lELd.NSAHPl 




r 
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EEPEODUCmiTY OP TOE 

OBKSNAL PAGE IS POOR 


LP-_1 NM AT( I Xil iiP. > GO TO 1 20 

NM » NI^AK 1) 

nn nn 

SAMPU-lflEL) * SAPP(I,IEU 

U2-E1I-1,1£U = SUEIUIEU 

no OIFF( I-l ,IEL> » DIFFd.IEU 

1 pn K MA H 1 -1 > « KMAT n > 

NMAT(5) « C 

1- Cr.MlMUE 

PRIKT in^ 

NMATCH s 0 

nn IP f «; 

fJAKEA * (N*l »*IAREA1 

n » ?. n » sukT 141 sq?pS I 

ISIZE = A5IZF (N,n 

I D » 2»1M(D/2 .0* 1. Cl - 1 

PRINT lOA, NSIZE(N). NAREA. Dt ISIZE, 10 

10 riMA Tr H » N»^ATCH ■»■ NSIZEIM 

PRINT lOA, NMATCH 

BP-IJLRAI 


c 

SUBROUTINE MATCH ( ST .NRl . NHA T . S AMP .S 1 ZE . D If P. T HR ESH . I EL . J S > 

C 

C COMPUTE CONTRAST FOR 5 SIZES OF DISK OVER A 5 X 5 PIXEL AREA 

C 

INTEGER ST(H. SAUPIS.lt. SIZE(S.1>. SUMlll). DISK. BOKO. DSAVE. 

.BSAVE, CDISK, HSN(25,10>, VSN(25tlO>, HSD(25.10)> VSD(23.10) 

D NMATll I . QlFPiS.ll. KiSIIOi. KAREAtlll. ASI2g 

.(5,2), MSEJ(SOO) 

cnMHr.n /5E2HCP/ N<>.HS»j,vsK. HsnTV^n.!<;PL).KagF t.A^i7F,MS,M«;pQ 

c 

-C COMPUTE SLM.S HP nATfl tfflIIIF S OVER ANMUHR RIMGS 

C 

JSl « JS - 4»NR1 - ? 

KEL « KAREA(ll) 

K « 1 

SUM(l) * 0 

0Q. J Q5-JS£Q«1.KE( 

IF ( JSEG.NE.KARE A(K)»1 I CC TO 305 

K « K ♦ 1 

SUM(K» » o 

305 ^LH(K) » SLM(K> » 1 AB S < S T ( J SI ♦ 1 S P Q ( JS EQ I > I 

C 

-C LPlIP over _5 512ES OF DISK AND S X S PIXEL AREA 

C 

AO ^ C.r 

00 18 KSIZc»l,5 

JS « -IS - 

NO » N<(NSIZE) 

kB » N SIUilZ EtS I 

C 

on 1 R KHO Hrl.S 

JS * JS ♦ NRl 

nn 18 KVER»1,5 

JSl * JS ♦ KVER - 3 

LE_IK. «FR.Nr.l 1 GP TO ?0n 

IE (KHCR.NE.H GC TO 2 01 

X 

C COMPUTE DISK AND BORDER SUMS AT FIRST POSITION AND SAVE FOR RECURSIVE 
X CALC U LATIPN S , 

D.J-S1L_=_0 

BORO » 0 

OQ .MST7F 

DISK « DISK * SUM(N) 

21 an &Q— «_aoa.a-i.. S UM N SizEtNi 

BORO * BJRD ♦ SUM(2*NSIZE»1) 

DSAVE » ni SK ^ 

BSAVE » BCTO 

GO TP. <;?in 

C 

X CD.MPUTF RErilftStV E SLM ALONG ROWS 

C 

201 D DIS K « P 

UD 3A JS£Q*1,N0 

2A DDISK « DQISK ♦ 1 Afl S ( S T( J SI »HS N( JS E Q .NS I Z E H 1 - 1 AH S ( STt J SI »HS D 

.1 |5F.j,).MZF||l 
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35 BPRO « BHKO ♦ I A B S ( ST( J SI -»H SN USE Q , NS I ZE >5 M 1 - I ABS ( ST ( J SI *HS 0 

.( JSEU,kSIZc>5» H 
f 

C SUBTRACT SUM QF TERMS ADDED TO PREVIOUS DISK SUM 
r 

BORO * BORC - DOISK 
0 S A VF - 01 SK 

BSAVE « 0ORD 

r^n jn Qjifk - 

C 

r roMPiiTP cpriip<;T\/P m ai nwr. rni 

C 

^nn nnfCK * r 

DO 32 JSEQsltNO 

nni<;K s nnr^K ♦ i a a ^Tf j «> u f 7 f i n - i arsi ♦v^n 

• ( JSEQf KSIZEU ) 

DISK * DISK ♦ I/DISK 

DQ 33 JSE«J’l,r<B 

31 ^ncp a luTRI ♦ I& PS( 5T( SI J*^EQtN5I 1 1 - lARSISH JSt*WSD 

.( JSEUfNSI ZE^EI » 1 
RHRn r RHRr - 00 I SR 

c 

r F|WO QP 0|FFfDE^Jrp PfTMEEW AU fB Af. F n I <1 K TFRM ANR AVEBAf-F 

C BORJEK TERM 
f 

9210 CONTINUE 

(JIF = DISK/ASl 7E INSIZE. 11 - BLJRD/ AS I ZE I NS IZ E . 2 1 

IF (OIF.LT.ACI GO TO 18 

An a p 1 F 

Kl « KVER 
KJ * KHOa 

K3 = JSl 
KA s NSI2E 

18 CONTINLr 
r 

C STORF MATCH PARAMETERS IF CONTRAST IS GREATER THAN THRESHOLD 
C 

IF (AO.LT. THRESH ) RETURN 
NMATI Kl 1 I Nl ATI Kl 1 ♦ 1 

NM = NMATIKll 

SA MPI Kl .NM ) = I F 1 - 5 ♦ K? 

SIZEIKI.NMI * KA 

UIFflKl.M'l » An 

C 

C IF «MCH FOUND. SET CENTRAL REGION TC NEGATIVE DATA VALUES FOR 

C SKIPP ING ALGG.^ITHM 
C 

OC 10 JSE0»1,MS 
JJ » K3 ♦ MSLuIJSEJ) 

10 STIJJ) • I SIGNISTUJI ,-l» 
RETURN 

END 


• 
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c 

SUBROUTINE SCi<TSL (A l. A2 . 11. JJ) 

C 

C SORT flRRflV ai AMD ARRANGE ARRAY A2 CORRESPONDINGLY 
C 00 FROM ELEMENTS II TO JJ 

X 

DIHENSICN AKKt A2(1U IU(16lt IL(16t 

INTFCFR A l . t?. . T? . TTl. T T2 

LOGICAL*! SL 

SL » .TRUE. 

CO TO 1 

X — 

ENTRY SORTLS ( Al » A2 i II, JJI 

C 

SL » .FALSE. 

1 NO « JJ - II * 1 

M*1 

L*_U1 

JsJJ 

5 IP(I.GE.J) GOTO 7Q : 

10 K*I 

1J«I JtIJ/2 

T1 * AKIJJ 

T2 » A2(1J» : 

IF (Aim. LE. Tl) GO TO 20 

AKIJ) » Aim 

A2(IJ» » A2(H 

Aim » Tl 

A2(n * T2 

Tl » Ain J 1 

T2 * A2( I J) 

?Q L=J 

IF (Al ( JI.GE.TH GO TO AO 

AI1_UI._^1UJ 

A2( Ij» = A2( J> 

Al (J) » T l 

A2(J» * T2 

T1 = A 1 ( t J 1 

T2 = A2(l J» 

LL-IAlXlJ.lF.Tll GO T O AO 

A 1 ( I J » » Aim 

A.2U JJ_-- A211J 

Aim » Tl 

fi?m *-i 2 

Tl » A 1( I J» 

IL2-« -A 2 ( I XI 

GOTO AC 

30 AKLI = AKKI 

A2(L» » A2(KI 

AllKl X TTl 

A2(KI * TT2 

AQ L«L-1 

IF (Al (LI.GT.TII GO TO AO 

TTl » AKL) 

TT2 • A2(L» 

«;n K»K»i 

IF (i1 IKI IT Til •' 1 m eo 
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I U ( M ) = I 
I*K 

■■■■■ 


GOTO 8C 
IllKl s JL, 


IU(M)* J 

-J-=J 


^Qin - 80 - 


70 M=M-1 

1L1M-.AE.01 GH TC 7-5- 


IF (SL) RETURN 
JJD.2.---MD/2 


11*11 
1 2 « J J 


on no i*i»ND2 

-XL— r-^XlXlJ 

T2 * A2( I 1 ) 



A2(lll - A2CI2) 
A 1 ( T ? 1 s T1 



A2( 12) » T2 

M = I 1 ♦ 1 

no 

12 * 12 - 1 
RF TLKfi 
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I* 1L(M) 
J» !U( M ) 



80 

1F( J-l .GE. I I ) GOTO 10 
IF ( I .EC. I n Giro 5 




1*1-1 



on 

I » !♦! 



IF(I.tC.J) &CTQ 70 


T2 * A2( 1*1) 

IF 1 A 1 ( 1 1. 1 F. n t (.n Tn 

90 



K=l 



100 

AKKfI ) : AlIKl 



A2(K+1 ) = 42(K) 

K = K-1 


IF m.LT.AUK) ) GO TO 
AlIKFllaTl 

100 


A2(K^1 ) * 12 

GOTO 90 
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APPENDIX B 


Listings of the following data handling subroutines are given: 


SETREE 

CELIB 

CLTREE 

RETREE 


selects training data from the feature vector file 

determines the frequencies and feature vectors in the 
foun-dimensional histogram analysis 

classifies trees according to the majority of pixels 

reconstructs the tree templates at the match points 
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n’r> o 


^ SUt<-< lUT jrT.<f!E ( X , S, * . N-tCCS t N’5A *’P. K * t :^* ^') .. 

SFLECT_ G-r,Ui.'.; 1-vUTr) TlLcS .. _ 

oiMCf.s:i;N s(.'.'it'.SA.-.p) i. 

I tiTe r. '-■K sm .-.f . <i Z:' s 1 1 

L-'.GlCi.L TaPLT( ^5»25,5> 


iro fup«aV ( 10 15 ) 

ni_FljWAI_( 12 0. V GCCJPI,rO_ CF.LLS'..n_0. !_.rGTAL V) 

FOKMAT (U,2i;it>» 

10b FriF.MAT ( El 6) _! 


00 20 h=l,5 


a.?. = Et.O «■ C. + IJ / 3.1Gl=92t5 
1 A 3 £A » 0 


DO 10 IM.Nfi 
_Al_r. -L r- 


DL) 10 

aj = J - ( IU»1 >/2 


TiiPLTl I = .false. 

IF ( a 1 »» 2 * A 01 .K 2 > CC TO 1 0 
I 4 i<EA = lAivCtt *• 1 
J ’’If Lrn.J .u) = . TRUE. 


'lO Cni.flMjE 
^0 PI^1'<T 1C2. 1A.3EA 


J:’:'.. *i^ 51 LL 


fc'Etri (Tc» s 
l^O W^lTE (90 ' UL > S 


on 260 NC»1.2 
NTGT(NC) = 0 

NT(NO_f_0 

KS' =' C 
JC( 1 I » 0 


C '^XTPACT CATA F'i-< A TkEE 
f 


250 Rj_An (5. ir :.t.<'J=220.)_ L1/-.E . SALP. SIZE. 
■■ NT(' 7 ct * ♦ 1 


Lll'^' LltiE - 

b..,— AF ’.Tim .~'T ~ 

a E A J ( 9 T L n ) s . - 

SI 1 * SA"P - jK/2 
Dj Aa NS 1= l.’lx 

”’lF ( . t.OT . T 'PLI (NLl .i;Sl »S IZF ) ) 0"' TH AA 

N I = N 1 ♦ 1 


f:'i 6C ’nfT^i.a 

63_X( I L- j« ’(FJ ’ - 

4 4 s'ii '= 'si i ♦ 1 

AO Lll = Lll ♦ 1 _ : . 

NTOT(NC) = '4rOT(NC) ♦ M 


WCIT^ (A.1C6) '.tC.G). ‘-'C . LIVE. SA"P. SIZE 


C o'eTER'm inf" cell' NUMttEKS 

c 
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— — 

— 

~ — — — — • - 


* c' 

— 

SoB^'ioTPJc CFLlS (X, ‘St •-'Ct JC , KSt L':j) 

. . . 

. c 


DT«eK'snM x<Nu,i), t-»c(nt jc(i) 




LOuICALRl oYTc(A) 

cfJUIVALPNCf (,CJ. u'YT-(in •. 


c 





1^0 

D.J lAO NF= UN J 

riYTFClFJ = X(NE.NSl) _ 


c 

c 

CHECK FfiS CELL JCCUPAWCY 


c 


on 36 K=1,KS 



36 

IF (tJCU.Ej.JCI K) ) t,U TO 37 
COHTINUE 


c 

c 

START NEW CELL 


c 


KS = KS f 1 


“ JC(KS) = NCO 

HC(KS) » 0 ^ — 

c 


K * K S 



37 

MCIK) = MC(K) ♦ 1 
PETJR^I 



END 
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r 


SllPR(iUTlNE CLTREc IS. ST. NR. lAt-EAl. prCc , NApCS. 

•IS AMP) 


c 

c 

CLASSIFY DETtCTEO TREES ACv-CROING TP THE E:AJi •ITY OF PIXELS 


c 

c 


IF M.;Oc*l, READ CCU^D I NA TE S ^RJM JAPE UNIT R 
IF .1J,7E«2, READ CCuR Jl.'.ATi S F^Ol CARDS 



r 


OIMENSirjN ISEJI5.3CL). NCLASS(i). MCL AS S 1 E 1 . PCTISD), KSIE.3) 




INTEGER STm. SInR.I), PIXLISI, SAMP(EP), SIZFIEP) 
Di’UBLE PRECISIJN CLASSK3) 

, CL^15SC3) 




DATA CLASS /‘DC. 'Ck', 'SC'/, CLASSl /• DCCLINE'. 
. • sniL'/ 

' HEALTHY', 



100 

FORMAT ( *1 *.23X, 'TREE CL AS 51 F I C AT I DN ' / 2 1 X , 19 1 ' * ' ) // 5X , ' L INE NUMBER 
.' ,5X, 'SAMPLE , SIZE. CLASS, PERCENT'//) 



101 

102 

FORMAT 1315) 

FORMAT ( '1 ' , 30X, 'CLASSIF ICATICN SUMM AR Y ' /3 IX , 2 2 ( ' • ' 

)//20X,' CLASS', 

,8X, *S IZE NUMBER' ,8 X , ' S A MPL ES ' 1 10 X , * P FP C E NT '/20X , 5 ( ' • • ) , 9X , 1 1 ( ' * ' ) t 
.8X,7( '*'),10X,7( •♦•M 



103 

104 

PJRHAT I/A26.1 n,l 17.F18.2) 

FrRMAT (/AiX, 'TOTAL' il 1"/A1X,15( '*') ) 




105 

106 

FCRMAT ( 112. 7X.6I 16. 13 .A6.FA. 1 ) ) 
format I 15X.6I I6.13.A6.FA. in 

7 


c 

107 

FORMAT (52X, 51 '-') ,11X,7( '-')/! 56.F1 9.2/5-2X.5I '-') , 

11X,7('-')//) 


c 

c 

DETERMINE FIVE TEMPLATE SEQUENCES 

. 




00 310 NT»1,5 

R2 * INT>1 )*IAREA1 / 3.1AI59265 

. 

NS60 * 0 
on 20 1=1. NR 

II = I - (Nk^n/2 
on 20 J*1.NR 

JJ*J-(NR*n/2 

IF ( n**=2*JJ*"-2.0T.R2) CO TO 20 

NSEO « NSEu ♦ 1 

ISFO( NTtNSEQ » » NR*(J-1» ♦ I 



20 

^10 

CONTINUE 

PI XL I NT) * NSEC 



c 

c 

INITIALIZE ARRAYS. READ INITIAL CDURPINATES 



c 


PKINI 100 


on 1C L»l,5 
- Bi: .10 N f jj. 3 


10 

Ksa.in * 0 

or .22.. Jy6 = i 



_ 

»EAC (10) (S(J-.ItL)f IE L = 1.NSA?-P) 

..L - - - - 

* 



IF IMOOE.tu.l) REAO 18) PfEC, NT, 
JF ("nE.E..2) REAU (5.1PU HR?C, 

(SAMP(I) , SIZEII), I«1,NT) 
SAmPJI ) , SIZE (1 ) 



MIDRO* * ,,R/2 ♦ 1 

TtR.E C 2 NR EC 5 . - „t*iL M 1 0 P 0 N _r . . 1 



c 

Ui: 130 MPEC*illDRC«,NREC2. . .. . .... 



_ 

IF (NREC.NF.MREC ) CFi TO 95 
-DQ 7? iir«liilT ^ 




c 

L ja_EACH CLA5i_. 

C 



c 

SUiiRIJUTINL KCTRtt (S, 5T . r.R, lARFAl, NtrCl, fSAtP) 

c 

c 

HFCONSTRUCT TFt TKEt TEKPLATE AT THF HATCH POINTS .. 

c 

INT’^GER, S( NR.n. STl 1) , I Si 0 ( E , TOO ) P 1 XL ( F J , i AKPI5C ) , SlIFI 50 » , 

,CLA5S(5C) 

LOGICAL MATCHS . . - - - - - 

c 

c 

DETERMINE FIvE TE1PLATF SECUFNCES _ ^ _ 

c 

DO 2 NT* I, 5 

R2 » (NT^n^IAAEAl / 3.1A1592f5 

NSEQ » 0 

DO 1 1*1, NR 

A1 s I - (i,ia+l)/2 

DO 1 J»1,NR 

AJ » J - (NR*n/2 

IF ( A I**2* AJ**2.GT.R2) GO TO 1 

NSEO = NS£d ♦ 1 . 

1SE0( r,T,\SEO) * NR*(J-l) ♦ I 

1 continue 

C_ 

2 PIXLINTI » NSEO . ’ ' 


C IFAD IFMTIAL CGOROINATEif SET BACKGROUND TO CLASS '3* 

i 

READ (8» MRtC. NT, (SAMP(I), SUE(I), CLASSII), 1*1, NTI 

MATCH S * .TRU E. 

MIDRDW * NR/2 ♦ 1 

NREC2 * WKFCS - Nk » MIDROW - 1 


DO 7 J*1,NR 
90 7 1EL*1,NSAHP 


c 

1 

SIJ.IEU » 3 


c 

r 

WRITE OUTPUT FILE OF RECONSTRUCTED TREES 


DO 10 NREC »MI0R0W,NREC2 
IF (NRFC.NF.HREC » GO TO 95 

DO 180 MI* 1,nT 

lAREA * P IXL I S IZ F( HT) ) -- 

JS » NR*(ShMP( MT )- l-NR/2 ) 

DO 7F NSFjsl.I AR = A_ _ ... 


75 

JJ » JS ♦ lS=j(SIZE( HT) ,NSEO) 
STIJJI » TlASSCHTI 


c 

130 

CONTINUE 




IF (HATCHS) READ (8, END* 30) HKEC, ^T , 
. I s 1 , NT) 

SIZE! I), CL ASS (I), 


30 

GJ Tu 95 

MATCHS * .FALSE* 


C 

5 Hll A .1 W S ^ - fyJk . - . 



6 

Du 6IFL»1 ,NSAHP 

SIJm-I.IEL) * S(Jil.flEL) . 



A 

on ft IEL*1,N5AMP 
S 1 N R . 1 F 1 1 X ft 


r 

10 

WRITE (U) (S(1,IEL), IEL = l,NSAHP) 

• • 



DO 200 JW*lt»4i< 

•i^RlTE (Hi ( 5iJdl 1 1 EL ) # i£L*l i.N S AMP ) 

• 

V 

✓ 

KETU3N 

E14 0 

... — - . . 

I 

t 
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c^’^RODUCBILITY OF THE 
npimNAL PAGE IS POOR 


APPENDIX C 


The following subroutines were used in performing the supervised classi- 
fications: 


SUBLPC 
PTEST 
PC LASS 
MAUCA 
EFFECT 

SNOPAL 

NTEST 

NCLASS 

NOPACA 


Compute Gaussian statistics of training samples 
Test maximum likelihood classification of training samples 
Handle I/O for classification of the data 
Maximum likelihood classifier 

Determine class to be tested for linear classifier using 
training samples 

Compute discriminant function coefficients for one class 
Test linear classification of training samples 
Handle I/O for classification of the data 
Nonparametric (linear) classifier 
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c 

SUfi KUUTfNE StWUC (LC. HC^IjS. CLASS. 3. gH. fK. ft. NW. Mill 


1 - j u r c:rs ij cu i n . ■«! .Mb rftKfinricKa 

c 

DlKfNSICN rtvdl, M5(r.f'.3), EM(K'N.HM). EK(NN.NN. HMI. BfNN) 

DOUBLE P?tC15lCr) CLASSO:^), S(KM,'JN), DET 
LCGICAL^l LCC.'j.l) 

4553 FaR;>.AT ( *1 */2jX, •£ STlrt>TfO GSUSSIAN PARANE TERS */20X ,29 ( ' *' ) //5X , 
. . , .'MEAN VPCT..KS* ,1CX. *C.:VA UAMCE MATRICES') 

45S4 FuK.MAT (=15.2.( 5x, 16E7.?) ) 

.... . ._ .^5'7 FtiRPAT (/20X, 'O'Tt LHI.JBNT »'.1PE10.3) 

45 77 FORMAT ( // 1 2 0 . A U , I 6 , ' SAMPLES') 
f 

C COMPUTE MEAN VECTORS AND COVARIANCP MATRICES 

c . 

RRITf (6,4550) 

.. . Orj 11 10 13*1, MM . 

NCll • hS( Ii,t) . - 

NCIP » NS( 

OG 93C0 11*1, .d ) 
=M(Il.Ii) » L. 3 

DC 4000 \'5i*r<c ii.ucic - 

. 4<'C3 fM(ll,I3) = CM(I1,I3) ♦ MC(f;Sl)*LC(Il.NSn 

9C.ro FMdl.m = F.Mdi.Ij) / NS(I3,1) 

. .. . c 

00 9200 11*1, .,:c 
D.- 9>00 n*i , I 1 

'R( 11 .d.U) » 0.0 ~ " 

DC 9130 NS1 = I,CII..<C12 

91'’0 FK(11,I2,1*) » ER(Il,I’,m ♦ MC(MSl) * ( L C ( 1 1 , NS 1 ) -E M (1 1 , f 3) ) 
; « (LCi iz.hsn-FN( u.i i) ) 

EK( 11,12.13) * EK(I1,I2,I3) / n.S(IJ,l)-l) 

S( 11. 12) « F^( n.t 2, 13) 

S( 12, m * S( 11,12) 

92<'0 CONTINUE 

WRITE (6.4577) 13. CLASSIlJ). NS(13.1) 

30 5003 11*1. NN ” • " 

: SE'CS PJsJ.rjT.. 4554, E1(II,I3), ( FK ( 1 1 , 1 2 . 1 3 ) . 12 * 1 . 11 ) 

C INVERT.. Cj]yARlAI,C£...><AT.RICE5, CMPUTE GAUSSIAN FUNCTICN CONSTANT TERMS 

• ^ . 

CALL GASI..V (i. JET) 

on 1000 Ii*l,t,N 
QC iC^CQ J,:*l ,:#N _ „ . 


ItCO cK(ll,I2,I3) » S(I1,12) 


tiL13±^_r3«.S-»_iNNyftLCS( 2,.C*3^141 59265I ♦ ALOCt SKCL ( OCT 111 

lllf' tiHITE (6, <.557) OET 

B£TUi<') 

END 


C 

SUb>’i:iJTlNF PT£ST (IC « . \S , CLi^ SS. ‘=K, y, N N. 

-C_ .Cl.»SS.lI.IJES..Kr4H>N. OflTA SA^P^fS - P AKA '.TTC I C CL A5S ! F IC AT IO!J 

C . 

piH'NSILN PC ( 1 ) t XI iZ), '’(1), 

'^LUBL? PPcCISIl!. CLA5S(.‘.p) 

LJ(^JCA L->1 LC( r. ,.1) ^ 

iCA FCkm,sT (//iox, IsHAifiixAuL iCCUr<ttCY * tFfc . 1 , » h PFXC5 M/ 3 ''X , 3^ ( IH* I) 

?C_I8_FCRJL®.I.-.< ’ft <:.?.( ' * • ) / 3; X, •_» P'" SULTi pc '“Laj.c jp JC4T I Cri »• /30X t '• 

. TRAIMhG SAiPLcS * ' / 'X , ’c ( i <. i f/ // j , • NU''B F R CF VJUX9ER 

! ^RCENJ _SA'--PLC S_CLA SSIF IcD A S ' /6X , • C L Ai S SAM 

.PLFS Cn.<iitcT CJ^RCCT* ,1 .AV/U3X.1 CA9> ) 

?f^09 Ft.RHAT (/K,A9.I7,1 1C,F1?. 1 ,1 C I 9 , ( /A2X , 101 9) ) 

C 


PRINT 2003, CLA4S 

TF » C.O 

Oi l 1301 NC « l. hM 

OC 1221 KR=1,MM 

1221 KSlNhl » 0 

i\SC = NS(NCtl) 
?tCil = NS ( NC.<:I 
NC12 * NSlNCti) 


DJ 14C0 NS1«NC11,NC12 

DO 15^0 tlF=l,MN 

15C0 X(I„F) = LCCNF.MSl) 

CALL MALICA ( X . ICLA S St KS . EM, F K, P. , 1, NN, MM) 

1400 KSdCuASS) * KS(ICLASS) ♦ KCCNSll - 1 

£ 

EFF * 100.0 » KSn.C) / NSC 

PRINT 2009. NCt ClAS^INO, NSC. KS(NC). FFP, (KS(NCl), NCl=l.hM> 

I’Ol TF = TE ♦ EFF 

AJI F = TF / F L J AT(F!M 1 

PRINT 104, AVE 

R FTUftN 

ENO 
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•y-f>rO' 




c 

5UBRjUTI,NE PCLftSS C S .MC L 4 S S .Cl ASS . E M. EK,B. WREC S . NS AMP , NK ,MM ) 

C 

C r t aSSlFlES U . 1 KNJ..N D AT A _S AM PLES - PAAAMgTRIC CLASSIFICATION ; 

C 

DlMENSnN S (\:j.NSaMP >. MCLASS( NSA MPl . EH(1). EK(H. B(l>. KS(2Q) 

0GL8L5 PRECISION CLASSdMI 

2: iQ FORMAT ( *1 « /3 OX. 29 ( » » « 1 /?0 X . » » RESULTS OF CLASSIFICATION »«/3QX. 

. '•',7X,*DATA SAMPlES'.'»X,'*'/30X,29( **')///20X, 'CLASS ' ,15X, 'SAMPLE 

. S'.15X t_LP;pCiNr*/20X.E( '*!_1»2-U 9JL» 7( '»') ) ) 

2011 FORMAT I /I IS , A 9, 1 20, F2 1 . 2 1 

^012 FORMAT I /1<,X. 1 SMTCTul S A MPLE S . 1 20/ 14X . 1 3 1 1 H* > ) 

C 

on 1A15 KC*1,M1 ^ 

KSINC) « Q 

C 

: O.Q_J 2_N5.E£ J.NAJE C 5 

: CALL MALICA (S. MCLASS, KS . EM, EK. B. MSAMP. NN. MMl 

: U hRITE till MCLAS S 


C 

NTt lT « NSAM P * _NRJc,CS 


PRINT 2010 
..Dn._i501 i<C * 1,KM 

PCT * 100.0 • XS(NC) / NTQT 
IjOl PRINT 2011. NC, CLASSIMC). KSINCl. PCT 

PKIWT 20U. NFOT 

RETURN 
* END 

. • 



• • 


• 


, * 
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SUP«<GUT1K5 v .ftLlCt (XI. K^^A X^K S 


C fAXInOK tUcLIHb?D CL AS S I F 1 C ATI CK . . . _ 

C 

... OIl^FNSinN .XlCJS'fJ.iS) , KX4X(!>.SS)» , 6 K ( NN , NN , F O , b C' A ) , 

.0X(20). KS(l) 

DATm GMAXX /ZFFFrFFFF/ 

C 

. DCL_2.0aLV5JrlfNSS 


G>AAX > GHAXX 
on 19 00 I^l.iSW 


G - B(Ii 
OQ F.3CQ II = l ,NN 


OX( I 1 ) > Xl( 1 I .NSl ) - EK( I I . I ) 
SUM « 0.0 


00 6200 JJ^l .1 I 

A « DX( _tKUJjJ. Ull 

6200 SUM * SUM - A 

63CC C = G » (SUM » 0.5»A) » CX(ll) 


IF (G.LT.GMAX) GO TO 1900 
KMAXINS n = 1 


CMAX = G 
1900 c n N TI MUE 


2000 KSIKMAXINSI) ) 
BETU 9 '^ 


KS(KMAX(NSin ♦ 1 


END 


n- 

» 1 
• i 
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c 

SLP9CUTlNr EFFtCT t LC . XC , t.'i, CUSS, *'uC t DE t >-Nt M) 

C 

C cPFECTlVE FIGUkE Q F MgitlT F FATURF_S.E.L E(|T 1 Jf i* 1 T FP I 

f 

riyptiSi r’i ■■cm. fc( 2 ?) 

LCOlCAL^l LC(fJi^.i> 

JGuqLF PKcClSlj.. CLfliS(FK) : 

lOO FOOKAT ( ' 1 '/2aX, *£FFcC TI V= FIGUKFS O' ^FS I T •/ , 26 ('*')// ) 

14 3 FOFMA T ll S.A^y.'X , 1 6r 7. 4 / ( ! =X 1 1 6F7...4) > 

I'C FORMAT (/2CX, 'CQMJ INFD FIGIkEF CF MEK I T • /2 CX ,2 F ( • • • ) / ) 

IFl FORMAT ( 125. 49. F 1 G. S . 5 X . 1 3 1 F ) 

C 

c cthputes inter-class and intra-class distahcfs 

c 

print ICO 

IF INMl.NE.l) GO TO 20C0 

DC 1005 NF 1«1,NN 

on 2 I1«1.HH 

NCll « NS(11,2» 

KC12 * NS( 11.3) 

DO 2 12«ltll ^ 

NC21 > NS( 12.2) 

NC22 » NS( 12.3) ^ 

0E( II .I2.NF1) » 0.0 

00 3 LKl = NC11.NC12 

IF ( I 1. EQ. 12) NC22 * LKl - 1 ' 

DO 3 LK2 » NC21.NC22 ^ 

3 DF( II .I2.NFI ) « Oc(Il.I2.NFl) «• I ABS ( LCI NF 1 ,LK1 ) -LC ( NF 1 . LK2 ) ) 

.»MC(LKl)»y.ClLK2) 

2 OEI I2.I l.NFl) < 0CII1.12.NF1) 

1C35 CONTINUE 

DO 11C9 NC^l.MM 

11P9 HOC I NO » NC : 

C 

C C1MPUTES THE NORMALIZED FIGURE OF HFRIT OF ALL REMAINING PATTERN 

C CLASSES ALONG EACH OF THE FEATURE OIRECTIUNS 

C 

2fG0 00 loro J3*N»l.MM 

CFMI J3 ) » l.Q 

13 » MjC(J3) 

U3_S_^XiU-3jJJ 

00 3CC0 I*1.NN 

F C*^»N « l. J g, 

NST * 0 

SLIMl . C.O 

SUM2 » C.O 

c c"c:m^pute suMi’-~TriTAL of~int‘er Class’ distanTeI—frcm' class u t’oTall 

C '^E^AlMNG CLASSES 

c 

0j3. i)_ ,J A_« N « 1 , h M 

IF (JA.EO.JA) Go TO 6 

ii».. »_Mac( 

AU > NS( 14.1) 

NST » NST ♦ NSI14.1) ! ! 

SUMl * SUKl * OEdi.14.1) 

C COMPUTE SUM2 - TuTAL dr'DTsTA'N^l'TMON^ALL R*E^ CLASTES. 

• C Ea uiVAi F>| T TO 1N T:<ACLA SS DISTANCE UF ALL FF,\Aru ?!G .CLASS £S_.CO.N ^ LOS'^.E.g. 


c xs class _ . _ 

c 

Pg 7 

IF <J5.Fii,J3l TJ*;** 

If * MnC(J5I . .* 

su;^? * Su;-,2 ♦ )F ( K. i ) 

7 CwUIi^UF ^ 

c 

C Cg^PU T E KlNIMo.n rlG u^E uF :«r.RIT I NDIVIDUAL CLAS55S 13 AND 14 

Q 

SI = 0E(I3tU,l) / IAI3 ^ AKL-. 

^ S 2 * CE(Ii.l3,I) / ( AI j^(AI^-UO) ) ♦ DE(U,I4,I) / ( A 1 4* ( AI A- 1 . 0 ) ) 

F * 51 / 

'If’CF .LT.FCKI.4) FC.'il., = F 

4 CCNTTMIE 

C 

ANSI = r.ST _ . 

SuHl * Soil / (Ali^’AiST) 

SU>2 * Drilj.I^.I) / ( fli :i»(AI3-K g) ) ♦ SUM.2 / ( A K$T » ( ANST-1 ,0 ) ) 

~"FC(1) * FC*U.\ wJ.il / SU’^2“ 

rc( I ) » cAP( >i,o/rv:( I ) ) 

c 

_C CrrPUTP CF h I SP t ly JL^l lit and order by cfh to determine 

C THE KLST SEPmumoLE cLASS 

c 

icro CFH(J3) * CFM(jj) ’J' rC(I) 

CFWCJ3) * CF.S(jn (1.0/LN) 

^RITF '(6^l43) lit CLASS(I3). (FC(i^F), NF»1,NN) 

IC OO CONTINUE 

’call SURTLS (CFiit MuCt NWl , MM) 

.c 

WRITE (6.150) 

:_0Q_J251 NC ; 

NC »*>l*OC(.'Ki ) “’ 

ysi hJRITE (6.151) NC. CLASS(NC). CFM(NC_1J 

CALL* SQRT5L (M6C. riL’C , NKl^l. HH) 

_KETURN_ _ 

’ end* 


I 


T> 
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c 

SUBROUTINE SNDPAl ILC t HC »VS .CLASS t W .PCCtS tV t B.NWl ,NN , m M ,KN 1 tM MI ) 

C 

C SUPERVISED NON-PARAHETRIC LFARMK& 

C 


IHI^ENSION rtCm, . WCfM ,r.Nl) ,_fiC , Ydl, S(ll 

L':crCAL*r'LC(.,:j, 1) 

CnUPLE PRECISION CL>.SS(^*^)» S(V*a,r.Nl), 


“ LfGIfAL TEST 

DATA N1 /ICO/ 

^ 99 FOkMAT I/1P6EIS.4) 

no FCRHAT 1 18 , III tIfl.AX, 1P7E1A. 3/(31X,1P7E1A. 31) 

101 FURI1AT 1/113, AIC,10X,1P7l1^.3/(33X,1P7E1A. 311 

n? FORMAT ( // 22X,2?( 'O' ) /?2X, '* CLASS ' .13 ,A10 /2 2X , 2 2 (.'*')// ' ITE 

.RATIUN NO. '.SX.'ERKoKS'.IOX, 'LINEAR DISCRIMINANT COEFFICIENTS'/ 
.A22,' CTrlER'/» 

216 FORMAT ( ' i '/ 10 X , ' JROER ED C LA SSE S ' , 20X , ' £ LE KE NTS PF THE OISCRIMINAH 
.T VECTOR'/lOX, 151 PCX, ’=( '*')/! 

220 FORMAT I/7X, 'TOTAL ERRORS', 15) 
C 

C INITIALIZE W. Y, AND 6 ARRAYS 

C , 

NW » MOC(NWl) 
NSiil > NSINM.21 

NSm 2 * NSINM.i) 
• NW2 * NWl ♦ 1 

DELTA * Nh/ACO.O 
DC 1112 NFA»1,NN1 

1112 WINWl.NFA) s 0.0 
J « 0 

NST2 * 0 

OG 1110 NC 1»NX1.MM 

NC > KCC(NCl) 

NST2 * NST2 ♦ NSINC.l) 

NSll * NSINC.2) 

NS12 * NSI NC .3) 

DO 1110 NS1>NS11.NS12 
J * J ♦ 1 

YU) *-1.0 

1 1 10 BIJ ) * 1.0 - - - 

r 

C CCf'PUTf INVEnSt uF AC Tn A N5P C5 c ) A UH5RE *A* AUCKENTEC MATRIX uF SAMP^Ei . 

C 

130 1*1 .NN 


bC liC J*i »\M 

Sf I iJ| *_ 0.0 V 

cn 13 1 KCl=.sfcl 

\c *j'riC(uci» 

f.Sll = K'S(.<C,2) 

b S12 « N SI ,CjJLl — 

no 131 NSi=NSi l.NSU 

K1 *.LC(ltNSU_. 

IF (J.NE.Mai Kl * Kl*cC(J,NSl) 

1.J31 SI I,J)_* _S( If.JJ_t HC.(NSU*Kl 

130 SlJ.n « SII.J) 

S(tlNl.NNl) » NST^ — 

CALL GASINV (3. NN 1 # DET) ■ • 

C DO NI iTpRATfjNS UF THE HU-KASHYAP ALfflRI THH » 'UNL E SS ALL CCEFFIcTeNTS 
C CHANGE BY LESS ..TrtiUi -QeiJA ^^/A_ P.EkC E NT 
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A 

■ 


P.<1ivT' 1C2. N... CLASSIf.n) » (LASS(N/< ) 
Dr. KAf 



TEST = .TAUE. 

00 iici I* i,.j.>rr" 

« AlLktl.tl 

J » 0 

r: 3101 :js i*;>i-.i ffiSr.3 

j = j ♦ 1 

» S( l.t.Nl 1 

14P 

no i'-.c 

A2 * A2 ♦ SI I,K)<'LClK.>iSl) 

2101 

hINkI. 1) - *I.Vrtl,n «• MCINSl)*42CAeSIYIJ)| 
DC 20''0 VC l = ii.w2iMM 

IvC * MOC(NCl) 

NS 11 s NS(NC.2) 

NS12 = NStriC.i) 

OJ 2000 NS1 = .\S11 .NS12 

j » j ♦ 1 
A 2 = Sd.t.NlI - 

141 

00 lAl K-l.iVN 

A2 = A2 ♦ SI JtKl*LCIK,NSl) 

2"*00 

rtlfj’.n.l) * -IlJ^l.l) - ■iCINSl) «‘A2«APSIYIjn 
IF IABS(«IU«1. 1I--0) .OT. APSIDELTA^nOI) TEST * .FALSE. 

1 101 
c 

CONT P4UE 

C CaKPUTc NEn UISCK IMIWANT VALUES AL'D CLASSIFICATION EKRORS 
C 

NERRl • 0 
I * 0 

00 ICCA KS1«NjW1 «N5W2 

1 « ! ♦ 1 

IF (Yin. CT. 0.0) nm = fain ♦ 2 .o*yii) 
Yin > hlNoItNNlI 

1141 

DO 1141 NF2-UNN 

YU) * Yin ♦ kI Nl. 1 ,).F2) *LCI NF2.NS1 ) 

1 r 

IF I YI 1 I.LE.0.0) t.ERRl - NFFRl »■ NCINSII 
YI I) = Yin - dI I) 

r 

NFl.-;2 = 0 

... 

“oo' ‘ircF ut i. 

NC * f re INCl 1 . ... • . 

NS 11 > NSINC.2) 

NSl? » NSINC.3) 

DO 1 ( Lb NSl»NSil.NSl2 

I • 1 ♦ I . .. 

1 45 

"if ivid.ct.o.o) ain = rid ♦ 2 .o*yiii 

VII) = WINUI.NNI) 

bo 145 NF2»1.NN 

Yll) a YII) ♦ wl NWl ,i.F2) ALCIU'^2 ,t,51 ) 

" K P5 

IF IVII l.OT.O.O) a rjPRR.; ♦ 'C INS 1 ) 

VII) »-YII) - BID _ . 

C 

PRINT ICO. n.uEX. NFR*.'.!, NFRR?, I W I Lv'l .OF A ) ,_ NF A* 1 , N'N 1 I 

l'^40 

Tf iVe's't)’"gu tj nio 

CfiNTINUf . • _ . 

1010 

c 

NERk « NEr k 1 ♦ !JcRR2 

.?RINT_t2C. . - - - 

IP INhUNE.imi RETURN 
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(6t4l6) _ 

on 1^20 = ■ 

MC » MUC(N*) 

l^cO PAINT ICl. NC. CL^SSiriC). (W(Nii.NFA), NFA*1,MN1I 

PRINT ICl. M uCI.iM) ,.CLAiliHaC(hM)J| 

RETURN 

fNP 


SH&.<,:UT!\£ NTEST <LC. KC , NS. CLASS. W. t'CC, MN, 


.y._‘ 5 S Lf 1 _E S... <jl C.iN_DA TA iAfPLfS 
Of.ENSI 


H .'.C (1 I . NSt • M .3 ) . X(2C) 
LrCICtL-^l LCl^-M.l) 


M'NPARAI FTKIC CLASSIFICATION 

( 1 ) . K.. cm, KS( 2 r» 


DCUa LE PhECIsI - 4 C LA i S L‘:’? L 

(// JvX. iart4v£\A0c 


lOA FCiKHAT ( // ;>:x , 1 arti v£\AOc aCCJAACY *,FF.l,P|- o = ri C C M / 3 X . 32 ( 1 H * ) I 
2^08 F PFrAT I *1 */j 3 X, 29 ( RESULTS 0^ r(. 4 S S I ' IC AT I C N *'./3FIX.,'* 

. " ’T>Alf.lF;v,'SAMPLES P''/30X,29I •*'»///iAX,'N»)>DEK OF NUH8EK 

. PEKCc^i.T NUM3ER C= SAMPLES CLASSIFIED AS * /feX , » C L ASS ^SAN 

CU^^iFCT' ,1 ''A9/(<.3X,10A9) I 


.PLFS 


CGRAECT 



£ V J'f 

c 

rLFNnMi rtiiifri.* 

NM > NN » 1 



.IMI * Hr* - 1 
PRINT 2CGd. CLASS 

TE * O.C 

on 1301 mc«i.M!4 __ 


1109 

DO 11C9 I«1,HN 
KS( I » » C 





NSC > NSINC.l) 
NCll « NS(SC.2) 




c 

NCU * >iS<NC«i) 



DO lACC NS 1>NC11 ,NC12 
‘ Do 15C0 NF*l.NN 


1500 

XIN-FJ » LCINF.NSl) 

CALL NGPACA (X. ICLASSt .«> 

MOC, 1, NN, 

NNl , MMl » . . 


1400 

c 

KSilCLASSI * KSIICLASS) ♦ 

HCCNSl) 




EFF * 100.0 * KSINC) / NSC 
PRINT 2C09. NC. CLASSINC). 

NSC. K.SfNC) 

, EFF, (KSINCl), NC1»1.HM1 

# 

1301 

TE « TE ♦ EFF 

AVF = TE / FlOATHM) 



PRINT lOA, AXE 
RETURN 

FND 


• 


H' 
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'^^‘RODUcmiLrrY of this 

^viiNAL PAGE IB POoA 


c 

SUttBPCTlHE McLASS is, M CLASSt CL A SSt Mt MPCt NRECSt NSAMP, NN . HN) 

C 

A^S'FJ3_IC CLASSIFICATICN 

C CLASSIFIES DATA SAMPLES At.O il'CRFS CLASSIFICATION RESULTS CN TAPE 

C 

DIMENSION S(NN.NSAMP), MCL ASS( NSAHP) , U(l). MOC(l), KS(20) 

DPUBLF PREC l S ION CLASS (J1 MJ 

2010 FpR.'<AT ( *1 '/3wX.29( '*• I/30X, •* RESULTS OF CLASSIFICATICN *'/30X, 

. »» »t7X» '04TA JAMPL': 5 • f 3X , ' »»/?CX.2« M »»*) ///20X , 'CLASS * il5X, ‘SAMPLE 
.S' ,15X, •?EkCE;JT'/20X,S( • * • ) ,2 ( 1 5X ,7 ( • * • ))) 

__2 0 11 FO R MAT 1/1 1 7_, A 1,_1 2 0 1£ 2 2 . 2J 

2C12 FORMAT ( /i jX. I jHTOTaL SA MFLE S , 1 2 0/1 3X f 13 ( 1 H*1) 

C 


K M » N K ♦ 1 

ME 1 « eTf - 1 

J).r _1 C C9_ NC_* 1 . i M 

ir09 KSINCl -6 

£ 

PO 12 .\RtC»l.NRECS 

ajAp .1 lO LJ' 

I 

CA LL N OP AC A. ( S, . HC LAS S , . KS . W. MOC. NSANP, NN. NNl. MHl I 

12 WRITE I 111 MCLA SJ 

C 

_NJOT ■_N$ANP_*_NRECS 

PRINT 2010 

DD_J221 NC*1.H.M 

PCT » ICO.O • KS1.NCI / NTOT 

12 21 PR I NT 20 11 » NC j_CL A S 5 1 NC 1 1 KSINCl . PCT 

PRINT 2012, NTOT 


RETURN 


suesruTi NE kopaca (x, uv , ks# "t ’loc, nss, nn, mni» mmd 

c ■ 
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00 20 NS1= I.USS 

Dil_l„N«l = l 

G * In ( N S' I > N n 1 ) 

DO 2 NF1 = 1 

2 u - G ♦ n|-<Fliihl)^X(N'ltNSl) 

IF ( G *0 T o" TJ S 

1 CONTIMUC 

Nw(Msi) = hJtir.Mi+n 

GC TC 2C 

-■l-NW.I.US 1 I = hi,C ( Nn 1) 

20 KSINiilNSl) ) = KSINnd.Sl)') ♦ 1 

RFTURU 

END 
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